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ABSTRACT 


Presented are results of some studies to develop tools useful for the analysis of Venus 
surface shape and roughness. Actual work was focused on Maxwell Montes. The analyses 
employ data acquired by means of NASA's Magellan satellite. The work is primarily 
concerned with deriving measurements of the Venusian surface using Magellan stereo 
SAR. Roughness was considered by means of a theoretical analyses based on Digital 
Elevation Models (DEMs), on single Magellan radar images combined with radiometer 
data, and on the use of multiple overlapping Magellan radar images from Cycles I, II and 
III, again combined with collateral radiometer data. The goal is to develop tools and 
methods to support the generation of planetary surface roughness maps. Such maps can 
be used to infer information about the ages of volcanic units through the use of erosion, 
weathering and deposition process models. The generation of roughness information 
requires precise knowledge of the terrain slopes as well as calibrated images showing 
SAR backscatter clean of the effects of slope and dielectric constants. We focus here on 
methods for surface shape measurements. We discuss data issues and show numerical 
descriptions of surface shape. The subsequent derivation of roughness has only been 
defined theoretically. Implementation and assessment of accuracy will remain the subject 
of future work. 
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1. INTRODUCTION 


Surface roughness of planet Venus from NASA's Magellan Mission is an important 
element in the analysis of Venus geology [Arvidson et al., in press]. NASA's Magellan 
satellite has of course produced a tremendous wealth of roughly 4000 radar image strips 
about the planet's surface, but it also generated surface topographic measurements by 
means of radar altimetry, observed the surface's emmissivity through radiometer 
observations and performed sophisticated gravity measurements [Saunders, et al., 1992; 
Pettengill, et al., 1991]. This collection of data has provided scientists with information 
useful in modeling many of the geophysical and geologic processes on Venus; it has been 
speculated that this data exceeds the quality and completeness of similar data sets 
currently available about our planet Earth [Leberl et. al., 1991]. 

The roughness of the surface can be defined in terms of detail at a continental or global 
scale, or it can address the local small scale variations over areas in the meter or 
centimeter range. "Roughness" is a property of the terrain that is either described as a 
root mean square variation of elevations in a given area, or as a root mean square 
variation of slopes. One has the dimension of length, the other of angle. While the 
global or continental scale assessment of roughness may be based on altimeter and stereo- 
based topographic Digital Elevation Models (DEMs), it becomes an issue of examining 
each individual pixel to derive observations of local surface characteristics. However, 
this analysis requires that accurate high-resolution topographic models be developed in an 
area of interest, that SAR images be calibrated, and that radiometric information be 
included in the analysis. The use of stereo imagery offers significant support for this type 
of analysis. This leads to topographic shape information and enhances the interpretability 
and usefulness of the original data sets. Derived data products can be defined and will 
provide a critical input to the detailed study of surface characteristics and to conclusions 
about important geologic and geomorphologic parameters such as the age of lava flows. 

We will focus in the following report on the technologies currently available and 
applicable to the Magellan full-resolution image data, either in the form of individual 
products of each orbit, the so-called F-BIDRs, or in the form of derived mosaicked F- 
MIDRs. We will discuss the currently available and implemented methods to measure 
the surface shape and to extract from it the roughness of a surface on a global or 
continental scale as well as methods for sub-radar-wavelength measurements. 
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The combination of such DEMs with the source images leads to radiometrically modified 
images which have the effects of surface slope removed from the image's brightness data. 
We describe a method to accomplish such a product, with complete implementation, as 
yet, ongoing. 

The resulting DEM can be used to describe, in terms of slopes, the "roughness" of the 
surface modeled in the DEM. However, this is not of great interest. Significant is the 
question of the local surface roughness within the area covered by a single image pixel of 
size 75 m x 75 m. Methods to address this local roughness include the consideration of 
dielectric properties of the surface within a pixel, and attempt to exclude the effects of 
this property on the image pixel's brightness. We will discuss the theoretical model under 
which this dielectric property can be modeled and removed from an image. The resulting 
product will represent the radar reflections as affected by surface roughness. 

1.1 Approach Based on Stereopsis 

In [Tyler and Simpson 1992], [Senske et al. 1992] and elsewhere it has been pointed out 
that large scale slope data are available from Magellan altimetry. However, high- 
resolution terrain information is only available through stereo techniques while repeat 
track interferometry is a promising method for surface shape reconstruction, but is not 
supported by the Magellan data set. As [Farr et al. 1993] report, about 21% of the surface 
of Venus has been imaged in stereo, offering an opportunity to not only generate high 
precision topographic maps but to also create geometrically and, more importantly, 
radiometrically corrected image products. Images from which the radar reflection effects 
due to terrain slope have been removed are called here "calibrated". Such calibrated 
images provide back scatter data that are essentially a function only of the small scale 
surface roughness and the dielectric constant of the surface’s material. As [Arvidson et 
al., in print] and others have pointed out, Venus is a planet with considerably intriguing 
dielectric properties of surface materials. Stereo-derived topo data can support analysis 
of these properties. 

The focus of our effort has so far been to produce DEMs of high accuracy and resolution 
from stereo images, depending on the limitations imposed by Magellan's imaging 
geometry and quality. These DEMs are used to correct image geometric and radiometric 
distortions due to terrain, following the model proposed by [Coriander 1991], [Rignot, 
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1992] or [Bayer et al. 1991]. Included in these corrections should be normalizations for 
SAR processor errors due to original Pioneer Venus and antenna pattern projections as 
well as the use of the specific model of Muhlman's incidence angle dependence. 

The DEM results in a "calibrated" image which can be the basis of surface roughness 
estimates at each pixel. Separately, the DEM itself is a source for computing revised 
surface roughnesses. We will discuss in Section 5 the case of the DEM improving the 
SAR images, and in Section 6 the estimation of roughness parameters. 

We differentiate between two Magellan stereo-processing schemes. One is based on 
mosaicked images (F-MIDRs) and employs many approximations; the concern exists that 
the approximations invalidate the usefulness of the resulting DEMs for all but the most 
localized studies. The other scheme is based on the original range-Doppler observations 
of the SAR system and promises superior accuracy; however, the methodology and 
processing system must still need to be put in place and verified. 

As previous studies have shown, SAR images can be employed for so-called "model 
inversion"; the resulting surface power spectra help to date lava flows. On Earth such 
data were used by [Evans et al. 1992], [van Zyl et al. 1991] and [Farr 1992] to establish 
surface roughness. Transferring these methods to Venus will require a determination of 
detailed topography and an appropriate theoretical model. 

1.2 Approach Based on Surface Radiometry 


Magellan radiometer data are to be used to infer dielectric constants for individual lava 
flow units. A simplified view of the relationship between (Fresnel) reflectivity p and 
dielectric constant e is given by [Farr et. al. 1993] as: 



0 ) 


where we can then obtain an approximate value of emissivity as e = 1 - p, where e is the 
(average) emissivity. This dependence is in turn removed from image pixels to leave 
only roughness-dependent scattering in the SAR image data. These reduced images can 
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We also have refined the DEMs from F-MIDRs using various implementations of so- 
called shape-from-shading. However, in an effort to differentiate surface units by means 
of different surface backscatter properties, shape-from-shading needs to be employed 
very judiciously. Its use presupposes knowledge of surface backscatter, a factor we seek 
to determine. Therefore we speculate that shape-from-shading can only then be used for 
this task if we manage to employ three overlapping images with strongly disparate 
incidence angles. We will therefore need to deal with all three Magellan Cycles I, II and 
III in one joint analysis effort to resolve the ambiguities between surface slope, dielectric 
constant and surface roughness. Shape-from-shading is thus a process element for 
associating with each pixel a slope estimate, and integrating these estimates into an 
elevation model, while other process elements resolve the issue of dielectric and 
roughness properties. 

At the current time we have not implemented a system to employ three images. We also 
still need to add radiometry data to the image analysis. 

2. TOPOGRAPHIC SHAPE RECONSTRUCTION METHODS 

2.1 Terrain Models for Planetary Studies 

In past planetary missions the major source of terrain elevation data was altimetry. This 
is in itself a method of elevation measurement which is robust only when the surface is 
flat, such as with ocean surfaces on Earth. As soon as topographic relief becomes 
accentuated, altimetric echoes will be ambiguous and the resulting DEM will be 
inaccurate. This will be particularly distinct in the description of relatively small features 
such as volcanoes, craters, steep mountain ridges etc. [Leberl et al. 1991] have found 
examples of altimetric observations over such features that are in error by ± 1 km or 
more. In addition altimetry observations are typically feasible at large intervals of some 
kilometers only, leading to a DEM with widely spaced postings. While Magellan 
produced a Global Topographic Data Record (GTDR) with a spacing of ~5 km, this is not 
the true spacing of independent observations. These are in the range of perhaps 13 km 
and more, and support the interpolation of a regular grid of topographic elevations. 

Magellan's image coverage supports the use of stereopsis to extract DEMs at accuracies 
of ± 100 m, and thus will permit one to describe the surface elevations by a grid with a 
mesh size of perhaps 500 m [see Leberl et. al., 1992], This can be used for detailed 
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analyses of surface characteristics such as small steep slopes, stability of the surface, etc. 
As [Arvidson et al., in press] have argued, such detailed knowledge of the surface shape 
also permits one to unravel the surprising brightness variations in radar images which 
may be caused by temperature changes on the surface as a function of terrain elevation. 
These analyses need accuracies of relative elevation differences in the range of 100 to 
300 meters. 

A separate issue is the measurement of slopes over rather short vertical ranges of, say, 

300 meters. These have been the subject of work by [Connors, in preparation]. Many of 
the slopes are close to the incidence angle of the radar system and get therefore fully 
compressed into 1 or several pixels. As a result, a stereo matching process will not work 
to measure such slopes. We will discuss alternate methods as used by [Connors, in 
preparation] or [Leberl et al. 1991]. 

In summary we have a number of data sources for extracting topographic information of 
Venus at various scales and locations. Some products already exist and can be directly 
used (such as altimetry) while others must be derived using stereopsis or radarclinometry 
(shape-from-shading) as described below. 

2.2 Stereo-Methods for Creating Topographic Data of Venus 

The stereo algorithms used in this study include those employed in conventional mapping 
as well as those found in 3rd-party image processing software packages. We utilize 3-D 
techniques of stereo reconstruction as well as simplified geometric methods. 


Two packages were used for elevation and DEM-processing, a system to process full 
resolution image strips per orbit, F-BIDRs, developed at JPL (Figure 1) [Hensley and 
Shaffer, 1994] and the Magellan Stereo Toolkit (Figure 2) [Vexcel Corp., 1994], for 
approximate computations based initially on mosaicked products (F-MIDRs). 

We thus use two techniques to reconstruct terrain, one based on F-BIDR data [JPL, 1990] 
and one based on F-MIDR data [JPL, 1991b], The first technique uses a detailed 
ephemeris and range-Doppler stereo intersection [Hensley and Shaffer, 1994], [Leberl, 
1990], [Dowman et. al., 1992]. The second employs a nominal ephemeris and simplified 
parallax computations to build the height models [Leberl et. al, 1992]. It would be 
possible to incorporate detailed imaging geometry and accurate navigational data into the 
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simplified parallax computations. However, the relationship between mosaicked image 
products and the radar system's detailed bursts of radiation is a detour that can be avoided 
when starting from F-BIDRs. 

Different methods have been proposed in the literature to make elevation measurements 
using overlapping stereo radar images. Authors include [Welch 1992], [Schanda 1985], 
[Leberl, 1990] and others. For the present study we compare two methods of differing 
complexity: 

least squares range-Doppler stereo intersection with detailed ephemeris and 
navigational corrections 

versus 

approximate parallax computation using a nominal ephemeris and no 
navigation corrections. 

Since the process based on F-BIDRs has previously not been described in any detail 
except for [Hensley and Shaffer, 1994], we will present additional details of the BIDR 
processes. In contrast, we have previously discussed the simplified process for parallax 
measurements with mosaicked F-MIDR data and therefore refer to that literature [Leberl, 
1990], [Leberl et. al„ 1992], [Vexcel Corp., 1994], 

2.3 Using Overlapping F-BIDRs 

We present stereo-processing software developed and utilized in this effort for F-BIDRs. 
While the F-BIDR-based process is separate from the Magellan Stereo Toolkit MST, it is 
compatible with it so that data can transferred into MST for visualization. Currently, the 
entire system consists of 4 parts: (a) reading of F-BIDRs; (b) image matching; (c) range- 
Doppler computing for the stereo intersection; (d) DEM gridding and image geocoding; a 
5th part is being developed for radiometric correction of the ortho-image mosaics. 

(a) Reading of F-BIDRs 
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This first software component has a largely administrative function and will not be 
further discussed. Details of reading and processing F-BIDRs can be found in documents 
of the Magellan mission [JPL, 1990]. 

(b) Image Matching 


Image matching is the process that takes raw input images, both in the form of BIDRs or 
MIDRs, and produces correspondence image points for a feature found in two images. 

Due to the nature of radar imaging the stereo partners are rather dissimilar [Domik, 

1984]; therefore it is useful that special methods be used to create useful overlapping 
stereo images. For matching by machine and for visual viewing one needs to constrain 
the imaging geometry to cases where the radiometric differences in the overlap areas are 
not excessive. [Domik, 1984] has described typical cases and specific imaging 
geometries that can be employed for optimal matching and fusing. The stereo imaging 
geometry adapted for Magellan was the result of an experiment during Cycle II of the 
mission [Leberl et al., 1992], This suggested that the incidence angles or Desired Look 
Angle Profile (DLAP) for the second look at the surface be inside that of the initial look 
in Cycle I. This was in principle implemented, with the exception of the orbits passing 
over Maxwell Montes. There the initial look angles were very small since the feature is 
at a high geographic latitude. Therefore it made sense and was feasible to have the stereo 
mate be taken at a shallower look angle over Maxwell Montes. Detailed look angle 
profiles for these orbits were described in [Farr et. al., 1993]. 

The fundamental trade-off that must be made when matching Magellan images is that of 
speed vs. accuracy. While automated methods have been shown to compare with manual 
methods to within ± 2 pixels [Leberl et al., 1994], there remains a need to visually inspect 
and accept the automated matches. This is due to the possibility of outliers or blunders. 
Manual inspection of dense automated measurements presents in turn a new problem in 
that the manipulation and viewing of individual points from a set of perhaps millions 
becomes prohibitively difficult. 

Image matching is a well-documented problem domain with extensive literature. We 
currently use two different matching algorithms. The Automatcher by [Frankot et. al. 
1994 ] is a hierarchical matcher of square mixed arrays computing a correlation maximum 
and a measure of confidence for each match. The match locations are regular in one of 
the overlaying images. 
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The second method by INRIA is feature-based, which searches for interesting image 
areas where strong matches are likely to occur. As a result, this method is robust where 
featureless areas may exist in the object. 

The match method of Frankot-Hensley-Shaffer [Frankot et. al., 1994] attempts to monitor 
results dynamically using adaptive "robust" weighted least squares (RWLS) estimation 
that we believe eliminates most of the blunders leaving errors that are primarily at a noise 
level comparable to that which an experienced stereo operator will encounter. 

A trade-off also exists between resolution of the feature to be captured by the matching 
process and matching errors introduced by noise in the images. After extensive 
evaluation we have concluded that matching can be performed optimally, for Magellan 
images, at a pixel spacing between 8 and 4 pixels, limited inherently by the signal 
strength of the multi-look image frames. 

The measurement accuracy of stereoscopic manual collection is typically viewed as being 
one of the critical advantages of manual collection in that it offers the ability of an 
operator to "fuse” a 3-D surface and place a measuring mark as opposed to 2-D matching 
or automated methods which are relying on variations in the image gray values . While it 
is likely that the automated matching performance on specific data sets can be tuned to 
improve upon or even match the performance of manual collection, this remains a case by 
case process highly dependent on the scene and imaging parameters. 

A significant limitation of stereo matching presents itself in highly compressed 
(foreshortened) slopes facing the antenna. Such slopes will become so narrow in one 
image that there is no texture that would permit a match with the second image to occur. 
However, it is these steep slopes that are often of specific interest and methods must be 
found to single out areas of such slopes which then could be reconstructed by machine 
using concepts such as those proposed by Connors (in preparation). 

(c) Stereo Solution 

This step takes measured image coordinate pairs in the F-BIDRs and converts them into 
measurements of slant range and Doppler frequency. Knowing the satellites position and 
velocity leads to the computation of a circle in 3-D space in each image. The two circles 
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produce the terrain point at their intersection. The 3-D cartesian coordinates need then to 
be transformed into a sinusoidal projection. Each patch of overlapping F-BIDRs 
produces a "cloud" of points. The stereo equations for nonzero Doppler [Hensley and 
Shaffer, 1994], [Leber 1, 1990] are 

r i = ( T - P i,T - Pi) 

f, J T - P - V >> (2) 

r 2 = ( T - P 2 >T - P 2 ) 

_ ( T - P 2 , V 2) 

^r 2 


since R; = T - P and r; = range; = |Ri| = V(f - Pi, T - P;). 

Here T is the unknown position vector of the target scatterer with 3 elements, fi is the 
Doppler frequency, r] is the (observed) range to the target for each of the two orbits, Pi is 
the spacecraft position vector for each orbit, R; is the pointing vector from sensor to 
target and X is the radar wavelength; i = 1 ,2. We have 4 equations to solve for the 3 
unknown elements of T. 


Here again is a trade-off: the more rigorous F-BIDR-method requires greater data storage 
and processing requirements, and navigational corrections that may not be available. But 
the stereo solutions offer the most accurate possible results in terms of absolute and 
relative planimetric and height errors, while the approximate F-MIDR-solution remains 
robust, simple, easy-to-use and available across the large amounts of mosaic data already 
issued to the Magellan science community. However, its accuracy may be low 
considering the noticeably large positional errors in F-BIDRS which will lead to false 
stereo-parallaxes and elevation errors. 

(d) DEM G ridding and Image Geometric Resampling 

The point clouds per patch of overlapping F-BIDR images serve as the basis for the DEM 
which is obtained by interpolating the computed points onto a regular map projection 
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grid. The user controls the output grid spacing and coordinate window in producing the 
DEMs. 

DEM gridding forms the resulting DEM by interpolating terrain elevations at regular, 
user- specified coordinate system spacings. The choice of grid-spacing is dependent on 
the match spacing and as a rule-of-thumb in mapping should be no greater than 3 times 
that of matching. For Magellan data of Venus however, due to the high-frequency terrain 
features found in many areas, we believe that this rule will not apply as strongly. In 
particular we have many highly foreshortened features. This should lead one to match as 
densely as possible, causing us to push the limits of matching to the 4-8 pixel level. 

The gridding algorithm is based on the Akima method, which is widely available through 
the algorithm collections published by ACM [Akima, 1978], This enforces continuity of 
first derivatives and integrability, important properties for roughness determination, given 
the requirement to work with slope and not elevation data. An additional property of use 
is that its interpolation model passes through each of the input points. 

This process augments previously used Magellan gridding schemes by (i) managing the 
computation using small blocks of data and (ii) performing image resampling during the 
DEM gridding stage. Given that the locations of DEM points are known both in the input 
images as well as in the output XY-coordinate system, it was seen to be an appropriate 
place to actually interpolate the image gray values for each input image. This process 
therefore produces a DEM in a user-specified output projection along with co-registered 
and geometrically calibrated SAR images. 

These geometrically corrected images will also have to be radiometrically corrected for 
the effect of the now-known slope of the terrain. 

2.4 Using Overlapping F-MIDRs 

F-MIDR products were designed to relieve the scientist from working with the complex 
and storage intensive F-BIDR products. Thus the thin orbital strips of 300 x 220 pixels 
are replaced with mosaics from multiple orbit swaths. The overhead of manipulating 
stereo models of entire planet-long strips in order to process aggregate areas on Venus 
using F-BIDRs is traded-off against a reduced accuracy and completeness of the mosaic 
products, however covering 8192 x 7168 pixels. 
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This stereo-process consists also of several steps as described in the following. 

(a) Matching 

Essentially the same algorithm described for F-BIDRs, modified to work on simple image 
pairs, is used to obtain matches for the F-MIDR pairs. 

(b) Stereo Parallax 

The match points produce parallax differences by assuming that the coordinates of the 
corner-points, as given, are correct (from so-called "dead-reckoning", using the nominal 
ephemeris) and thus allowing an absolute datum reference given that all images are 
projected onto the Pioneer Venus topo model. The emphasis is not to compute absolute 
elevations of the terrain, but rather to routinely determine the relative elevation 
differences with respect to the de facto datum point (which may be in error due to 
ephemeris uncertainty). The DEM should undergo a “calibration step where it is leveled 
to other data (at least 3 points) such as altimetry or BIDR stereo. 

The parallax differences are then simply converted into elevation differences by an 
equation which assumes known radar incidence angles, and it replaces the spherical wave 
front by its tangent plane. This approach is based on a number of simplifying 
approximations. Its DEM is that which would be produced from the more rigorous 
approach would the ephemeris be error free, at a constant orbit altitude at a constant look 
angle across the entire mosaic and without the spherical figure of the planet. 

The stereo parallaxes represent the surface elevation above Pioneer Venus topomodel. 

To obtain the elevations above a spherical datum surface it is necessary to back in the 
topomodel. 

As part of the mission's preparation for stereo data collection, an analysis of F-BEDR 
radar images evaluated parallaxes and concluded that cross track navigation errors can 
occur that cause up to 1 km error in the computation of height differences between points 
located in different F-BIDRs [Hensley, unpublished JPL-Mission material]. 

(c) Gridding and Resampling 
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A variant of Akima’s algorithm is built into the image processing subsystem IDL 
(Research Systems Inc.) and is being used as part of MST. It is being employed to 
generate MST-based DEMs and to geocode or rectify the SAR images. 

This F-MIDR related software differs from the F-BIDR gridding element, although both 
use Akima. MST's does not allow gridding of very large point sets encountered in typical 
Magellan stereo processing. The MST version is dependent on physical memory and has 
problems when processing point sets in excess of 50,000. This limit does not exist in the 
F-BIDR processing chains. 

2.5 Using Altimetry 

The direct measurement of terrain elevation is of course through altimetry. MST supports 
access to the global altimetric data from Magellan. This is the recommended data set for a 
broad overview of terrain elevations. Its validity is limited to areas of gradual topographic 
relief. In areas of topographic relief (where stereopsis is applicable), altimetry can 
produce erroneous elevation data, distorting particularly elevation profiles over features 
with delicate shapes such as volcanoes, craters and such. 

Raw altimetry profiles as observed by the satellite and not yet processed through 
interpolation to create the global data set may contain information about small features 
such as craters. As [Leberl. et. al. 1991] have shown, this can, in combination with the 
radar images, lead to improved measurements of the elevation differences in craters and 
volcanoes. Access to raw altimetry echoes is through a separate software element 
denoted MGMDQE and available through the MIT Center for Space Research. 

2.6 Shape from Shading 


A promising technique for achieving maximum terrain model resolution is that which 
employs radarclinometry or "shape-from-shading" to interpolate stereo-derived height 
measurements using the input SAR images and the slopes obtained from the stereo data, 
MST contains a method, described by (Thomas et. al., 1991 ], that is an expansion of a 
method previously described by [Frankot and Chellapa, 1987]. It is currently 
implemented for a nominal ephemeris and has been tested extensively on MIDR 
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Magellan data. A number of improvements to this algorithm, however, have been 
identified. 

The basic idea behind shape-from-shading is the hypothesis that an image grey value 
(digital number DN) is greatly dependent upon the incidence angle, and that therefore the 
incidence angle can be estimated from DN-values. Various authors have proposed 
algorithms to compute incidence angles, terrain slopes, and terrain elevations [Wildey, 
1978; Kirk, 1987]. These do not typically employ multiple radar images, which can be 
easily merged in stereoscopic analysis of multiple images. 

Shape from Shading (SfS) depends on an estimate of the surface backscatter properties, 
so that image brightness variations can be correctly interpreted as a result of slope 
variations of the terrain. The weakness of SfS is its inability to determine backscatter 
properties. Instead one simply must assume the backscatter characteristics in an often 
arbitrary manner and as a uniform item across the entire imaged area. The model used is 
either a cosine-reflectivity law or the reflection models proposed by [Muhleman, 1964] or 
[Hagfors, 1964], each with global parameterization chosen somewhat arbitrarily. 

Major improvements may include not just single Muhleman's or Hagfors' backscatter 
constants but multiple constants that are input as masks into the imagery and cause 
associated lookup table elements to be referenced during reflectance map calculation. 

The dielectric constants, derived from the radiometer data in a prior step, should be used 
in the above reflectance map computation. 

The capability of SfS to utilize precise navigation information found in the F-BIDR 
processing parameters and employ a precisely computed look angle off-nadir to replace 
the simplification of a constant look angle is needed as a further improvement. 

2.7 Other Methods 

Direct monoscopic or multiple image methods exist for making spot elevation 
measurements. [Leberl. et. al., 1991] reviewed the methods in existence. These may 
include the use of the observable length of shadows or layovers, etc. It may exploit the 
knowledge that a feature such as a volcano or crater is symmetric, and that therefore the 
two opposite sides of the feature can be treated as if they were the two images of a stereo 
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pair. Another approach is non-stereo exploitation of overlapping images, as [Dalke and 
McCoy, 1968] have demonstrated, and as Connors (in preparation) is now using to 
determine the slopes of small features which cannot be matched in stereo. 

3. DATA ISSUES 

Our work focuses on a small number of areas that are of interest for their geometric and 
radiometric qualities and because of the availability of special orbital correction data. We 
plan to compare various types of topographic datasets that cover the areas. At this point 
we have only performed work on Maxwell Montes. However, the other areas remain on 
the list of things to do as resources permit. 

3.1 Maxwell Montes, 65 North, 0 East 

This area (Figure 3) has had extensive previous analysis due to interest in both its steep 
terrain and surface properties [Pettengill and Ford 1991], [Alexandrov et. al., 1986]. For 
this reason Maxwell Montes has been paid particular attention by the Magellan Mission 
as well as the Magellan Project in providing special orbital coverage, custom data 
processing and navigation ephemeris corrections. For this region, specific corrections 
were made to the navigation data, provided by research on the north pole control network 
[Davies et. ah, 1992], and ephemeris refinement by [Chodas et. ah, 1993], and applied 
during reprocessing of the basic image data products (BIDRs). 

The select availability of these navigation corrections currently constrains the use of the 
most sophisticated and accurate methods of stereo solution to only Maxwell Montes. 

The specific site is presented in Figure 4, using a portion of an F-MIDR stereo-pair. 

3.2 Ovda Regio, 2 South, 74 East, 8 South 71 East 

This area is considered part of Aphrodite Terra. The area under study has been imaged in 
stereo under the Stereo Experiment using images from Cycles I and II. Given its 
proximity to the equator it affords an opportunity to compare differences 
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Figure 3. Maxwell Montes Cycle 1 F-MIDR covering Latitudes 63.5 to 67.5 degrees 
and Longitudes 357.0 to 2.0 degrees, encompassing 240 km x 450 km and showing 
the study site in the enclosed box. 



(b) 

Figure 4. F-MIDR Cycle I (a) and Cycle III (b) stereo image pair; coverage is from 
Latitude 65.3 to 66.2 and Longitude 356.9 to 0.75. Nominal spacecraft altitude was (a) 
1127 km and (b) 1115 km with look angles off nadir of about (a) 22 and (b) 33 degrees 




between the MIDR and BIDR results that is free from map projection distortions. 

This area has been matched manually, including terrain breakline acquisition, giving also 
the opportunity to assess the performance of the automatching (for mosaics) with respect 
to individual operator measurements for Magellan SAR images. Later in this report we 
will discuss the performance of this program as used on optical photography on Earth 
having known accuracies. 

Finally, this area has cross-polarized emissivity data, also acquired in an experimental 
mode. This dataset, as we shall see, can be uniquely used for roughness determination. 
This area is not discussed any further in this report. 

3.3 Artemis Chasm, 35 South, 145 East 

This area has been extensively evaluated by using stereopsis, altimetry and monoscopic 
methods as well as for geologic and structural content by Connors (in progress). For this 
reason we include this area as a study site for future work since the terrain complexity is 
great, offering multiple fault directions and contrasting high and low frequency terrain 
patterns. 

For now we need to defer study of this area until further navigation corrections and other 
ancillary and image data become available. 

4. SURFACE MEASUREMENT RESULTS 

A number of DEMs were derived predominantly over the Maxwell Montes test site. 
Before addressing roughness determination we describe the methods and parameters used 
to generate these height models and present a comparison of the results. 

We show results for the 3-D range-Doppler BIDR stereo solution of [Hensley and 
Shaffer, 1994] and that of the approximate parallax method [Leberl et. al., 1992], using a 
"nominal” ephemeris, for the same area and mosaic image data (MIDRs). Given the 
large distribution of MIDR products, it is additionally useful to assess the extent to which 
approximate methods produce results that differ from those obtained by intersection 
methods. We speculate that the accuracy of the approximate method depends on 
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geographic latitude, with errors increasing significantly with distance from the equator. 

For Maxwell the F-MIDRs have been reprocessed with ephemeris corrections so 
navigational errors should be minimized. In addition, for completeness and consistency 
we present Magellan and Pioneer Venus altimeter results. 

4.1 Stereo Measurements on F-MIDRs and F-BIDRs 

This section discusses the results for the test sites for which stereo processing was done, 
including use of stereo intersection, parallax and shape-from-shading, all applied to 
Maxwell Montes. 

A DEM was produced using the stereo processing software developed at JPL (Figure 5a) 
The images strips were matched at intervals of 4 pixels and the DEM was gridded to 1 
pixel [Hensley and Shaffer, 1994], The DEM was subsequently averaged to compare with 
2-pixel DEMs generated from parallax stereo. The averaging was done to decrease the 
noise in the DEMs. 1 

Figure 5b shows a DEM of the same area produced from the MIDRs (using about thirty- 
five thousand match points) using software that is part of the Magellan Stereo Toolkit 
[Vexcel Corp., 1994], It uses the nominal ephemeris and approximate parallax schemes 
described in [Leberl, et. al., 1992], Matching was at every 6 pixels while gridding was at 
4 pixels; the DEM was supersampled using bilinear interpolation to match the resolution 
of the DEM from F-BIDRs. It should be noted that the MIDR DEM is thus somewhat 
“smoothed” while the BIDR DEM is somewhat “noisy”. These are properties of primarily 
the match and grid spacings used. 

Both DEMs were generated using the automated matching algorithm described in 
[Frankot et. al., 1994], The matching software exists in two formats, one for 
conventional raster images (such as MIDRs) and one specialized for Magellan BIDRs. 
The accuracies and characteristics of the matching methods are discussed later in this 
report. 

Figure 6 shows a series of profiles through the detailed area for the parallax and 3-D 
intersection DEMs. The profiles are indicated in the images of Figure 12 by the black tick 
marks at the left edge of the Cycle III image. Wire-mesh perspectives are shown in 

1 This work was done partly at Vexcel Corporation under the current contract. 
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Figure 7 for the two DEMs and Figure 8 shows IHS images with color encoded as height. 
Figure 9 illustrates a perspective view of the DEMs with the Cycle 3 image overlaid. 
Figures 10-12 show more presentations of the results including SfS DEM contours over 
the Cycle III image. 

Table 1 lists statistics for the two DEMs. While large errors due to ephemeris have been 
corrected by the reprocessing of the MIDRs used in the parallax computations, the large 
discrepancy in the height values we believe is due to coordinate system projection errors 
as well as noise in the BIDR DEM. The mean difference between the MIDR and BIDR 
DEM after registration was 759m + - 457m. The minimum and maximum differences 
were -2288 and 2244 respectively. The MIDR-based version relies on a simple 
"counting” formula for image pixel location as opposed to the BIDR version which, for 
each image framelet pixel, goes through exact coordinate transformation during 
processing. 


Table 1. DEM Statistics for Maxwell Montes 



Resolution (m) 

min 

max 

mean 

stand, dev 

BIDR 

75 

4884 

10732 

8155 

1016 

MIDR 

300 

5272 

10081 

7747 

759 

MGAlt 

4.5k 

4687 

9563 

7614 

821 

PVAlt 

10k 

5270 

9445 

7437 

731 

SfsMIDR 

150 

5228 

10043 

7747 

759 



(b) 

Figure 5. DEMs from F-BIDR (a) and F-MIDR (b) stereo, presented as images 
with contour overlay at 250 meters. The area is about Lat 65.3,66.2, Long 356.9,0.8. 
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Figure 6.Elevation Profiles of the Maxwell Montes Test Site Using Different DEMs. 






Figure 7. Wire-mesh perspectives from the Southwest for (a) the F-BIDR 
DEM and (b) the F-MIDR DEM. Vertical exaggeration is 20. 
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(b) 

Figure 9. Perspective views of (a) F-BIDR DEM and (b) F-MIDR DEM. Vertical exaggeration is 2 with 
elevation angle 15 degrees, azimuth rotation -150 degrees. The Cycle 111 F-MIDR ortho from parallax 
stereo was used in draping. 



Figure 10a. Perspective of F-BIDR DEM as in Figure 9a but color encoded with the IHS image of Figure 8a. 



Figure 10b. Perspective of F-MIDR DEM as in Figure 9b. but using the IHS image of Figure 8b. 
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Figure lib. Simulation of F-MIDR DEM as in Figure 11a 




c . „ h of Shape- from-shading refined DEM from F-MIDRs. SFS iterations used 

Figure lie. Simulation as in Figures lla,b of Shape trom snau g 

was 10. 










4.2 Stereo versus Altimetry 


GTDR and PV Altimetry were obtained for this area (Figure 6). The values were 
extracted using a weighted Lagrangian interpolation that enhanced the smoothness of the 
altimetry at the resolution of the stereo maps. MST extracted the area but it seems as if it 
needed to be shifted by 225 150m pixels to register to the other DEMs. Table 1 shows 
the elevation ranges. 

4.3 Refining a Stereo DEM with Shape-from-Shading 

Shape-from-Shading (SfS) was applied to the images of Figure 4 using the MIDR DEM. 
Radar image simulations are shown in Figure 1 lc for one case. Figures 1 la,b illustrate a 
simulation of both the stereo MIDR and BIDR DEMs indicating an improvement in detail 
by means of SfS. Future work will try to exploit the ability of the increased resolution 
DEM of SFS to remove a greater amount of local incidence angle effects thus leaving a 
more precise backscatter dataset for use in inversion models. 

An elaborate set of results has been compiled from which Figure 1 lc is but a small 
sample. However, the analysis of the results requires further work. This analysis will be 
presented in the final report. 

4.4 Geometric Calibration of SAR Images using a DEM 

The SAR imagery is being corrected for the geometric effects of topography (terrain 
distortions). However, also radiometric effects need to be modeled and removed. Thus 
the accuracy of the DEM (from stereo or other means) affects the radiometric and 
geometric accuracy of the resampled ortho- image mosaics [Curlander, 1991], [Bayer et. 
al„ 1991], [Rignot et. al., 1992], 

During BIDR processing the stereo DEM is deformed by regridding the irregular 
elevation measurements onto a regular grid of postings through interpolation. Similarly a 
warping field may be derived by regularizing the computed shifts in image pixels due to 
the terrain relief that was determined in the stereo intersection process. Both MIDR and 
BIDR methods must produce ortho-image products that are registered to the DEM. 
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The above step is referred to as "ortho-rectification", "geometric calibration", or 
"geocoding”. Geocoding requires further that the corrected image be transformed onto a 
defined map projection. Sometimes the transformation of a geometrically distorted image 
into a map projection is referred to as “geocoding with terrain corrections . When this 
step is complete the location of pixels is assumed to be correct, while the backscatter 
values are still fundamentally based on incidence angles that assume a spherical planet 
shape normalized to the Pioneer Venus topo model. Thus local incidence angle must be 
accounted for in the correct computation of backscatter, a process that involves the 
simulation of the Magellan orbit as described in Section 5. 

4.5 An Experiment to Assess Automated Stereo Matching 

Using Non-Magellan Imagery 

In order to place the performance of the automated image matching algorithm used in this 
study contextually within classical image matching capabilities we have selected an aerial 
stereo photography pair over terrain near Albuquerque, N.M. We performed an 
experiment that compares the MST matcher with that of an experienced stereo operator 
using a standard photogrammetric stereoplotter to make manual measurements. Since 
photogrammetric methods of stereo measurements have been extensively tested and 
studied and our algorithm for matching Magellan images is robust in the sense that it can 
be made non-sensitive to large noise it seemed a worthwhile test to show the accuracy of 
the automatcher in an environment which provides stable and accurate measurement 
control such as that of conventional aerial stereo photogrammetry. 

Since the original manual measurements do not produce image coordinates, but only a 
DEM, we have to perform number of transformations from image pixel into fiducial 
space and then to actually generate a DEM from the automated match data. Results are 
reported in terms of differences between the DEM from manual collection versus that 
from automated. Thus errors in the transformation model as well as in the DEM height 
interpolation must be considered. 

We matched photographs represented each by 7K x 7K pixels at 32 pixel spacings where 
each pixel was 0.25 m on the ground as determined through the film scanning process. 
Due to the high variance of parallax in the stereo model we used a 6-level hierarchical 
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matching scheme that began with blocks of image data with 256x256 pixels and ended 
with 32x32 blocks. A total of -20,000 points were matched by machine. 

Results (Figure 13) give a mean difference of 12m and a standard deviation of .7 m. 
Expressing this in terms of pixels we find a random difference of about ± 3 pixels in 

elevation. Due to a base-to-height ratio of 0.6, an elevation error x is caused by a 
matching error ( x j. This leads to the conclusions that matching differs by +_4 pixels. 

This shows that there is some systematic bias where the matcher DEM consistently 
underestimated the terrain. The overall detail of the DEM and the ability to follow the 
terrain were very reasonable. Reasons for the bias may include photogrammetric 
orientation errors, match errors or interpolation errors. The difference image of Figure 
13c displays a trend of decreasing accuracy as slope increases in steep terrain and 
increasing accuracy at lower elevations. This is perhaps because the algorithm scans from 
left to right. An improvement might be achieved by averaging multiple direction scans to 
eliminate this artifact. 

It is believed that through parameter adjustments, increased point density, or improved 
scans, we could further improve these automatcher results to within +- 1 pixel. 
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Using Magellan Imagery 


We also performed a similar experiment using a portion of the Magellan Maxwell Montes 
test data. In this experiment we compared DEMs (a) produced within MST from match 
points using the Frankot-Hensley-Shaffer algorithm; (b) produced by an algorithm 
developed at INRIA (National Research Institute for Information and Automation) in 
France. The INRIA algorithm is feature based and utilizes the epipolar geometry for 
match point refinement. Figure 14 presents the resultant DEMs along with contour plots 
and the input SAR images. No quantitative comparisons were yet made but we note here 
that the correspondences were reasonable at 200m levels, with more errors apparent in the 
INRIA algorithm probably attributable to the sparse feature-driven point sets. Future 
research may address the possible merging of a feature based algorithm with cross 
correlation in a statistics based algorithm to perhaps obtain measurements closely 
resembling those obtained by hand, which can model terrain breaklines. 

5. REMOVING THE EFFECTS OF SLOPES FROM MAGELLAN IMAGES 

Removal of slope-effects in SAR images has been studied using a number of methods. In 
this report we focus on removal of local incidence angle area and antenna pattern effects 
for Magellan images (see for example [Cheng et. al., 19941) using stereo techniques. As 
an extension of this concept we are pursuing the interpolation of stereo elevation 
measurements using grey values (shape-from-shading) as a way to increase the resolution 
of the stereo topography and thus the accuracy of the roughness determination. After 
removal of processor area and antenna pattern local incidence angle effects we will be left 
with backscatter values that are functions primarily of small scale scattering and dielectric 
constant. Residual errors, according to [Rignot, 1992], consist of SAR processing 
calibration errors, as well as errors from the stereo and SfS processing. 
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5.1 An Approximate Computation of Backscatter from DN-Values 

To some degree we can simply view or segment the resulting reduced imagery, especially 
over uniform cover. For Venus, and Maxwell Montes in particular, we must further 
remove the effects of dielectric constant to get real roughness data. The degree to which 
incidence angle estimation errors affect determination of dielectric constant has been 
studied by [Rignot, 1992], He did show it to be significant for calibrated radar image 
data. 

For Magellan, we can get backscatter go, normalized to the Muhleman model from 
G() = 0.2 (pixel DN - 101) dB 

from which we can then correct the values directly using approximate ephemeris. A more 
accurate approach involves using the Magellan radiometric resampling coefficients which 
are in a range-Doppler coordinate system and then to actually form corrected DNs in a 
forward manner as described next. 

5.2 A Rigorous Algorithm to Compute Backscatter from DN-Values and Terrain 
Slope 

To reduce the geometrically corrected image mosaics we need to form the Magellan orbit 
and, using the high resolution stereo DEM, determine the look vector to each individual 
DEM scattering element. From the look vector and the local slopes we can derive the 
local incidence angles for each DEM posting (and thus each pixel element in the mosaic) 
which in turn determines the various radiometric corrections to apply to the mosaic grey 
values. In order to form the look angles the nearest bursts and their parameters must be 
located and accessed with respect to the pixels (scatterers) in the image mosaic. 

In order to apply radiometric corrections using actual incidence angles due to terrain we 
must first remove the radiometric compensation factor applied during processing of each 
burst. This factor, C rg (r,f), is a function of r, range coordinate, and f, Doppler frequency. 
C r g(r,f) is defined as: 

(47t)3N r 2N a R4(r,f) 

rg r ’ PT(t)KNpGBAQGA 2 (6e,0h)CrCdG r (t)?i 2 f(I)G s tate 
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where 


R(r,f) = slant range to pixel 

transmitted power of burst 

TW - transmitted power during calibration 

N r = range FFT length 

Na = azimuth FFT length 

Np = # pulses in burst 

K = transmitted power during calibration * sensor gain * processor gain 

GA(9e, Oh) = 2-way antenna pattern function, including loss through cable 
connecting sensor and antenna 

GbaQ = residual gain factor left out of the BAQ reconstruction process: 

n . 1-73 S 2 

g BAQ- 127.0 

where S2 is a function of the threshhold (a constant per burst) 

G state = gain for non-standard redundancy configurations 

0 e , Oh = elevation and horizontal angles off-boresight 
Cr = conversion factor from time to ground range 

p. - 

Cr -2sin(I m ) 

where I m is the MRP (mid-range point) incidence angle 
Cd = conversion factor from Doppler to along-track distance: 

Cd = Rm(- — 1 t> - 

2IV s lVl-(<Vs ,P >)2 

where Vs = spacecraft velocity, P = boresight pointing vector, and V$ and 
P are the corresponding unit vectors 

receiver gain of burst 

t->rv ) - rece iver gain during calibration 

X = radar wavelength 

I = incidence angle assuming a spherical surface 
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= 0.01 18cos(I) 

() (sin(I)+.l 1 lcos(I)) 3 

Note that the incidence angle I figures into C r g(r,f). In practice, a constant value Im is 
used per burst where I m is the incidence angle at MRP (the mid-range point) assuming a 
spherical surface. Furthermore, Crg is only computed at nine points per burst and a set of 
quadratic equations are derived to approximate C r g. These equations are of the form: 

C r g(r,f) = a(r) + b(r)f + c(r) 2 

where 

a(r) = ai + a 2 r + a 3 r 2 
b(r) = bi + b 2 r + b 3 r 2 
c(r) = c l + C 2 r + c3 r 2 - 

Thus to remove the radiometric compensation factors applied during processing, we must 
do the following: 

(1) Compute r (range) and f (Doppler frequency) of the selected pixel in the BIDR ortho 

(2) Compute C r g(r,f) using a(r) + b(r)f + c(r) f 2 
and 

a(r) = ai + a2 r + a3 r 2 
b(r) = b i + b 2 r + b 3 r 2 
c(r) = ci + C2 r + c3 r 2 

where the coefficients ai , a2, etc. are those used to radiometrically compensate the 
current burst. 

a n 

(3) Remove C r g(r,f) from a n (normalized backscatter) to get IX(r,f)l 2 = q ( r ”f) where 
IX(r,f)l is the pixel value after range and azimuth compression. 

Now we want to apply the correct value of Crg(r,f) to the pixel value I X (r ,f) I to get the 
new normalized backscatter On- 

(4) Compute the actual Crg(r,f) using 

(47t)3N r 2 N a R 4 (r,f) 

PT(t)KNpGBAQGA^(0e>6h)CrCdGr(t)X 2 f(I)Gstate 
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and the actual terrain (and thereby the actual incidence angle I). Highlighted terms are 
assumed to be constant over the burst. Details concerning this computation follow: 

(4a) To find 0 e , 0h, and I we need to find the spacecraft position X s (t), velocity V s (t), 
and boresight pointing vector P(t) 

(4b) Given X s , V s , P, compute Cd 

(4c) Given X s , r, and f, intersect the range-sphere and Doppler cone with the actual 
Venus surface model (iteratively) 

(4d) Once the intersection of the surface, range sphere, and Doppler cone is found, we 
can compute 0e and 0h and then find GA(0e> 0h) 

(4e) We must also compute I: 

(i) Form the vector from Xs to the surface intersection point Pt (position of target) 
and convert to a unit vector (see Figure 15). 

(ii) Find the slope of the surface in the direction parallel to a plane containing X s 
and Pt 

We have ^ and ^ of the surface at or near Pt using the DEM, where y is N/S 

dimension and x is E/W. The normal to the surface is found by taking the cross- 

. . de , de 

product ot ^ and ^ . 

(iii) Compute the normal to the plane containing X s and Pt 

(iv) Project the surface normal onto the normal found in (iii) 

(v) Subtract the projection from the surface normal 

(vi) Compute I as follows: 
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Figure 15. Definitions Used in Radiometric Corrections 

a) Negate the unit vector from the sensor so it faces away from it 

b) Compute dot product between the surface normal and unit vector to get I 
(vii) Now find f(I). 

(5) Now that we have GA(0e, 8h). Qb and f(I), compute C rg (r,f) using equation in (4) 
and apply to Gn = IX(r,f)|2 Crg(r,f). 

As a result, we now have computed a new image in which each pixel is represented with 
its backscatter value clear of the effects of terrain slope as derived from the DEM. 
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5.3 Implementation Status 


The program for radiometrically correcting geometrically calibrated mosaics for errors 
due to local incidence angle (and therefore slope) is currently being developed. This 
program represents the fifth element in Figure 1 which shows the F-BIDR processing 
parameter data and geometric mosaics as inputs and the radiometrically and 
geometrically calibrated image mosaics as output. Vexcel has implemented similar 
programs for the ERS-1, JPL Airsar and SIR-C SAR datasets and plans to utilize these 
code segments in finishing this program module. Other elements of this program include 
the Magellan orbit simulation as described above, and the access and manipulation of the 
processing parameter files, much of which is contained in the BIDR stereo processor. 

6. MEASURING ROUGHNESS USING A DEM 


A height model can be used directly to assess surface roughness at various scales. 
Statistics such as RMS height are only somewhat useful in that they do not provide 
information on the relative scales of roughness that may exist [Farr et. al., 1994], Instead 
investigators have adopted a number of theories for describing the scale and degree of 
roughness for natural terrain surfaces. 

6.1 Methods of Roughness Characterization 


In particular we are interested in models for diffuse scattering which is primarily what the 
Magellan SAR at off-nadir look angles of up to 40 degrees will represent. Diffuse models 
(Bragg scattering) assume that small scale roughness (smaller than the wavelength, ie. cm 
to m scale) scatters the incident electromagnetic energy in a non-specular manner. 


(a) Raleigh based roughness is measured by means of a root mean square change in 
surface elevation, Oh- 


As an indication of roughness for Magellan we can examine the Rayleigh roughness 
criteria (in centimeters) according to the Magellan wavelength of 12.6 cm: 


smooth : : 


X 1 

Oh < = cm 

25 cos9 2 cos0 
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rough: 


cm 


X 

ah > 

4.4 cos0 


2,8 

COS0 


where Oh is the root mean square variation in surface elevation. 


(b) Using correlation length, RMS height determination for a given scale is possible 
through methods as described in [Ulaby et. al., 1986]. 

(c) Alternatively spectral techniques that use a log-log plot of the power spectral density 
of a profile of elevation data file vs. the wave number [Evans et. al., 1992], [Brown, 
1987], [Farr, 1992] followed by a linear fit to determine slope and offset offer a way to 
generate data at different roughness scales. This method assumes that the power spectral 
density (PSD), denoted by G(k) of all natural terrain surfaces obey the empirical relation 


G(k) = C k ' a (4) 

according to [Brown, 1987] and others, where k = ^2L \ s the wave number and X is the 

distance along the elevation profile. This allows the extraction of slope as a measure of 
relative scale (steep is small cm scale roughness and flat is large) along with offset which 
has been described qualitatively as the “roughness amplitude”. The slope can be related to 
self-similar fractal dimension D, by D = 2.5 - a/2 if it can be shown that 2 < a < 3, 
offering additional methods for analysis and use of terrain data [ Brown, 1987]. 


An important element in this technique is related to the moments of the power spectrum. 
From Fourier analysis we find that the Oth moment mQ relates to roughness as follows 


( CjlKt 

G(o>) dco 

“0 


1/2 


(5) 


while the second moment is the variance of slopes and the fourth is the variance of 
curvatures [Brown and Scholz, 1985]. Detailed analysis of these moments leads to the 
conclusion that the fractal dimension changes with frequency and thus can be constant 
only over a limited range. 

For detailed information on small scale roughness however, some form of 
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microtopographic profile is required such as in [Farr, 1992]. It has been shown [van Zyl, 
1991] that the power spectra for such microtopographic models can be estimated using 
model inversion techniques that use calibrated backscatter values to compute roughness. 
The relation 

Oh = 40.0 e offset (6) 

was derived empirically by [Evans, 1992], where Oh is roughness (RMS height) and 
offset was derived from the power density. 

For this study, the DEMs alone are inherently limited in terms of the roughness 
information they provide, by their resolution, typically much greater than the radar 
wavelength of 12.6 cm for Magellan. In order to infer small scale roughness we must 
resort to the use of SAR pixel grey-value and other ancillary data as described below. 

(d) Roughness from slope values [Arvidson et. al., in press] defines "roughness" as an 
angular value using local slopes within a scattering area. This is denoted as CJ rm s ar >d 
represents a root mean square value of slope angle. 

6.2 Implementation Issues 

We plan to select one or more models for direct roughness information extraction for 
implementation under the Magellan Stereo Toolkit. Candidates include producing slope 
and offset images using elevation data and other readily accessible, generally applicable 
functions. The correct use of these functions must be the responsibility of the 
“knowledgeable” user. 

One exciting aspect of these roughness studies is the potential for placing stereo cameras 
or CCD arrays on specialized landers (eg. see [Arvidson, et. al., in press]) and acquiring 
microtopographic profiles directly of the venusian surface for comparison with the 
inversion models that are discussed in the next section. 

7. MEASURING ROUGHNESS USING A MAGEUUAN IMAGE AND 
RADIOMETRY DATA 

7.1 Methods of Dielectric Determination 
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For Magellan we have a direct means to estimate the dielectric constant by using the 
Magellan radiometer data. From the radiometer data, which measure emissivity, we can 
estimate the dielectric constant e from the formula (see [Arvidson. et. al., in press]): 


e a (0) 

where e a is the average emissivity which can be computed by averaging observed H and 
V emissivity values (for Magellan Stereo Experiment only); 0 is the local incidence angle 
and e is the unknown dielectric (assumed to be the real part only). Note that equation (7) 
simplifies to equation (1) if the incidence angle approaches 0. 


= 1 - 


1 cos 0 - V e - sin 2 0 

2 

1 | e cos 0 - i e - sin 2 0 j 

'cos 0 + V e - sin 2 0 ^ 


^ '£ cos 0 + V e - sin 2 0 


(7) 


Since the resolution of the radiometer data is much lower that the SAR we must assume 
continuity of the geologic units being measured; that is we will find aggregate pixels that 
are of the same chemical properties for areas that are larger than the size of the 
radiometer footprint. 


Another method has been reported by [Rignot, 1992], for JPL Airsar images, that is 
independent of “roughness” (using the Small Perturbation Model) 


Y = 


[(e - 1) sin 2 6+ e]| COS 6 + (e - sin 2 6 ) 1/2 / 
ecos (0) + (e - sin 2 0) 1/2 ] 2 


( 8 ) 


where y is the ratio between calibrated SAR VV and HH backscatter values, 0 is local 
incidence angle, and the dielectric constant is again to be solved for. 


We need to find a general way to compute a pixel-by-pixel estimate of dielectric 
constants for Magellan data (without cross-polarized emissivity data), which we expect to 
be quite high for instance for the Maxwell Montes [Ford, 1994], [Tyler et. al., 1992] area. 

Estimates of dielectric constant allow reduction of the grey-values to images that are then 
a function only of small scale roughness and calibration errors. 


7.2 Methods of Roughness Determination (Model Inversion) 
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(a) Using the results of ground based experiments without dielectric properties 


One way to convert from the qualitative reduced backscatter values to quantitative 
roughness indicators requires inversion of a scattering model to produce pixel-by-pixel 
maps of "offset" as described by (ie in equation (6)) [van Zyl, 1991], [Evans, 1992]. For 
instance, we can compute roughness, ah (RMS height), from an empirical formula by 
[Rignot, et. al., 1992], where a^HH is the backscatter of a scatterer at HH-polarization: 

a h = q - 2 . ^ hh . (9) 

(exp (-0.722) +o°hh ) (k sin (0; ) cos 4 (0i )) 

Here k is the wave number, 0i is the incidence angle. This formula is only valid, however, 
for calibrated SAR data and slightly rough surfaces but does offer a method dependent 
only on HH backscatter values. 

For Magellan there needs to be derived a form or extension for the inversion of the small 
model perturbation [Rice, 1951], [Kim and Rodriguez, 1992], [Ulaby, 1985] or other 
model that can be used to infer small scale roughness as described above. 

(b) Using a Model with Specular and Diffuse Components 

In [Arvidson, et. al., in press] we find the most detailed formulation to date on inversion 
for diffuse scattering of Magellan data. It is briefly stated below. 

First we form a specular component to the scattering model as we assume a two-scale 
(Kirchoff and Bragg scattering) linear combination approach [Kim and Rodriguez, 1992], 
We denote the quasi-specular component of backscattering as a^QS(0) where 0 is 

incidence angle: 

o°QS(0) = -- P - Q — N(0). 00) 

2 cos 0 

Here po is the Fresnel power reflection coefficient at nadir, calculated as 



and N is the facet number density (see [Arvidson, et al., in press]). Similarly, the diffuse 
component is modeled as 


Go,d = 2 po n cos n (() (12) 

where n is a controlling exponent for weighting the backscatter according to incidence 
angle. When we combine the two components linearly and expand N (see [Arvidson et. 
al., in press]) we get the forward direction of the inversion formula for roughness Orms as 
a function of backscatter oq: 


Co = ( 1 - f D ) Oo.qs + f dOo.d = ( 1 - f d) P ° f * — exp | tan^G 1 + 2 f D p 0 n cos n ()) 

2 tarn Crms ' 2tan z ORMSt 

(13) 

We have observed backscatter GO, fD is a percentage the user assumes to be diffuse; 
values for Orms are angular as defined in section 6.1(d). 

(c) Other Methods 

Other methods that are alternative to the SPM include UPM (Unified Perturbation 
Method [Kim and Rodriguez, 1992]) and use of solutions to the full wave equation, 
incorporating multiple reflections [Bahar, 1994], 

While external calibration does not exist for Magellan its internal calibration effects are 
fairly well known and it is a good candidate for the investigation and development of 
inversion techniques and correlation with surficial morphology for the determination of 
processes and rates for deposition such as that being investigated by [Arvidson et. al., in 
press]. 

7.3 Implementation Issues 

We plan to select one or more models of roughness determination for implementation 
under the Magellan Stereo Toolkit. Candidates include use of the two-scale inversion 
model of [Arvidson, et. al., in press], and other model inversions that are available for 
Magellan SAR data. A preliminary specification of modules is given in Figure 16 
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showing direct roughness determination, dielectric determination and backscatter 
inversion models as functions. We would like to host a suite of various methods that can 
be used and compared by investigators doing roughness studies. 

More theoretical work must be done on appropriate models and their inversions for 
Magellan data where HH and VV cross-polarized emissivity data is unavailable. Figure 
16 details the system specification for a roughness modeling process similar to that 
described in [Arvidson, et. al., in press]. 

8. MEASURING ROUGHNESS USING MULTIPLE MAGELLAN IMAGES 

AND RADIOMETRY DATA 

The most interesting aspect of using multiple Magellan images for roughness 
determination is the aspect of SAR imaging that inherently provides information about 
roughness of different scales according to the look angle geometry of the radar. The task 
remains to find ways to combine the information from the multiple images into a single 
quantitative description of roughness. 

Presumably the Cycle III data will have more variance after reduction [Bayer, 1991] then 
the Cycle I image since it has higher grazing angle but there may be areas from the 
steeper imaging that may be enhanced as well. Incorporation of the right-look Cycle II 
image would give an additional roughness image that should represent roughness at a 
slightly different scale than the stereo data. These represent differences in small scale 
roughness that are separately enhanced by the respective imaging geometries. 

9. OUTLOOK AND PLANS 

We plan to increase our test site coverage with more BIDR data and to finish implementation of 
the radiometric correction programs. This completes the 
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Figure 16. Data-flow diagram (DFD) for dielectric constant and roughness determination 
programs. Note: P1-P5 software programs are described in Figure 1. 


collection of an end-to-end Unix-based BIDR stereo processing program set that can be 
used by Magellan investigators in their detailed surface studies. We also intend to clarify 
and codify availability of ephemeris corrections and supplement them with relative 
corrections where unavailable. 


Also we will investigate, derive and if possible exercise one or more of the techniques for 
terrain roughness characterization including correlation, inversion and fractal modeling; 
in particular we will investigate an appropriate form for an inversion model for use with 
Magellan data. 
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Our estimates of roughness at various sub-wavelength scales can be input to process 
modeling programs that can infer rates of processes such as that done by [Arvidson, et. 
al., 1992], 


Finally, we plan to show that increasing the DEM resolution through a rigorous stereo- 
constrained application of shape-from-shading to F-BIDRs produces improved 
radiometric correction capability of the SAR images. These new calibrated images, we 
hypothesize, will in turn propagate into better estimates of surface roughness and hence 
surface process modeling and dating. 
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